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Abstract 



We show that the multipole vector decomposition, recently introduced 
by Copi et al., is a consequence of Sylvester's theorem, and corresponds to 
the Majcwell representation. Analyzing it in terms of harmonic polynomi- 
als, we show that this decomposition results in fact from the application 
of the harmonic projection operator and its inverse. We derive the coef- 
ficients of the usual harmonic decomposition from the multipole vectors. 
We answer to " an open question ", first raised by Copi et al., and re- 
ported by Katz and Weeks, by showing that the decomposition resulting 
from their corollary is unstable. We propose however a new decomposi- 
tion which is stable. We generalize these results to complex functions and 
polynomials. 

1 Introduction 

The recent CMB data have emphasized the necessity to handle data on 
the sphere, like the temperature [fluctuations] of the CMB on the last 
scattering surface. Since the role of Fourier transform on the sphere is 
played by the expansion in spherical harmonics, this approach is widely 
popular. The techniques are now standard and have been widely used for 
analysis and interpretation of the results. 

To analyze a function, the first step is generally a separation of scales 
provided by the multipole development / = /(^) . The sum is over all 
integers, but a cut-off is made at some L to take into account the resolution 
of the data. The spherical harmonics decomposition is completed by 



Recent analyses using various methods, have reported some 

signs of anisotropy and/or non gaussianity in the CMB data. On the 
other hand, such effects are expected in some theoretical models of the 
primordial universe. This motivates an active research for such effects in 
the present and future CMB data, for which the relevance of the spherical 
harmonics approach has been questioned. 
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In addition to the unfriendly behavior of the Yim under spatial rota- 
tions, this has motivated interest towards other techniques. In this regard, 
Copi, Huterer and Starkman (0, hereafter CHS) have recently shown the 
association of a series of £ " multipole vectors " to any multipole /(£). 
They suggest to represent /(^j by these vectors instead of by the usual 
firn- Having shown interesting properties of this set of vectors (in partic- 
ular the fact that they are rotated by the usual representation), they give 
convincing arguments for the relevance of their new techniques, whose in- 
terest is confirmed by application to the CMB data. This result was later 
confirmed and interpreted in terms of polynomials by Katz and Weeks 
hereafter KW). However, as mentioned by both groups, many ques- 
tions remain about the signification, relevance and interpretation of this 
decomposition. The present paper intends to contribute to this subject. 

In section 12.21 following a recent work by Dennis I show that 
the CHS result is in fact a direct consequence of the Maxwell multipole 
representation. 

In section 12. 3L I recall first some well known properties of harmonic 
polynomials. Then, I show that the correspondences established by CHS 
and KW identify to the harmonic projection and its inverse (whose ex- 
istence is established in this context). This allows a shorter proof of the 
decomposition of an harmonic vector, and provides the link with the usual 
spherical harmonics approach. In particular, this allows 12.411 to estimate 
the fern coefficients from the multipole vectors. This also allows 113. Ill , to 
give a (negative) answer to the " open question " of stability, raised by 
CHS, and reported by KW, about their multipole development; on the 
other hand, I propose a new development which remains stable. Finally 
||1J, extension to complex polynomials is discussed. 

2 The real vector decomposition 
2.1 Multipoles and harmonic polynomials 

Functions on the sphere may be seen as reductions to 5^ of functions 
on the embedding space IR^, in particular polynomials. 

There are many complementary ways to consider a multipole /(^j: 

• as an eigenfunction of [the Laplacian Ag2 on] the sphere with 
eigenvalues Xe = —I {£ + 1): I call f(e) an ^— eigenfunction; 

• as a function with definite squared angular momentum A^, since A32 
identifies with the squared angular momentum momentum opera- 
tor J^. A further classification, by the projection Jz of the angular 
momentum, leads to the usual Yim, as the normalized ^—eigenfunction, 
with eigenvalue m of J^. 

• as a vector of the (2£+ 1)— dimensional irreducible representation of 
the rotation group S0(3), or of its universal covering SU(2). 

• As the reduction to of an ^-harmonic polynomial: an homoge- 
neous polynomial of degree £ (hereafter an ^-homogeneous polyno- 
mial) of H^, which is IR'^-harmonic, i.e., verifies AP = 0, where A 
is the Laplacian on IR"^. 
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Hereafter, I note Harm(^) and Harm(^, IR) the vector spaces (on ]R 
and C respectively) of ^-harmonic polynomials with complex and 
real coefficients respectively. I note Hom(£) and Hom(£, IR) the vec- 
tor spaces (on IR and C) of homogeneous polynomials of degree £ 
(hereafter ^- homogeneous polynomials ), with complex and real co- 
efficients respectively. 

• As a completely symmetric traceless tensor of rank £. The natural 
development of any ^— homogeneous polynomial P £ Hom(i?), 

P(X)= Yl F,„,,,...,,,x'^x'\..x^\ (2) 

involves a completely symmetric tensor Fij^^i2,...,ii = J'(ii,i2,...,if)- 
Thus, any homogeneous polynomial correspond to such a ten- 
sor. For an harmonic polynomial, this tensor is traceless: any 
£— harmonic polynomial correspond to a completely symmetric trace- 
less tensor. The harmonic projection (see below) is obtained by 
taking the traceless part. 

• An equivalence class in B[om(£), see below (12.311 . 

2.2 Multipole vectors and Maxwell representa- 
tion 

Let us consider the /-homogeneous polynomials P of the form 

P{X)=A{X-ui)...{X-ui), (3) 

where ^ is a real constant and the Ui are unit [and real] vectors of IR^, 
i.e., points (^directions) on the unit sphere. I call 

Vec(/,IR) C Hom(/,IR) (4) 

the set of such polynomials (which are in general non harmonic). 
CHS have shown the one to one correspondence 

Harm(£,IR) >^ Vec(/,IR) (5) 
H 1-^ vecH, 

with vecH of the form (^J above. 

They explicit the correspondence in the tensorial notation; symmetriza- 
tion of the tensor Ui^ Ui^-.-Ui^ gives the tensorial form of vecH. Then, the 
traceless part gives the tensorial form of H . This allows to solve explicitly 
the correspondence (although with heavy calculations). 

Following 151, we show that this correspondence is a consequence of the 
Sylvester's theorem. (Note that the classical proof of Sylvester's theorem 
implies Bezout's theorem, like the proof in |S]. However, |2] provides a 
different proof). The latter states that any /-harmonic polynomial with 
real coefficients can be uniquely written as 

H{X)=r^'+^ V„iV„2-V„,-!-, Vi/ G Harm(/,1R), (6) 
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with = ^ X ■ X, the Ui as above and the directional derivatives 
= Ui ■ V . This is known as the Maxwell multipole representation, 
and [21 has shown that it implies the unique decomposition 

H = vecH + r^ Q;Q eRom{£-2,TR). (7) 

This establishes the correspondence namely the CHS's result. 
To resume, 

\/H G Harm(^,]R); 

HiX) = r^'+^ VmVu^.-.Vu, - ^ A {X ■ ui)...{X ■ ut) + r^ Q; (8) 

r 

Q G Hom(^- 2,]R). 

The proposition by CHS, and then by KW, to characterize any H 
(or, equivalently, /(^j) by the constant A and the 2£ components of the 
Ui (rather than the 2£ + 1 a^m's) offers the main advantage that the Ui 
rotate under the IR'^ representation (which corresponds to £ = 1). As 
shown explicitly by this corresponds to a factorization of the [spin] £ 
representation into a product of ^ = 1 representations. 

Now we show now that this decomposition corresponds to the har- 
monic projection. 

2.3 Harmonic projection 

The well known decomposition (see, e.g., 

Hom(^) = Harm(^) Hom(^ - 2) (9) 

means that any ^-homogenous polynomial P G Hom(i?) can be uniquely 
decomposed as 

P = HP + Q; UP € Harm(^), Q G Hom(£ - 2). (10) 

The £- harmonic HP is called the harmonic projection of P 0, and Q is 
{£ — 2)-homogenous. It can be checked that 11 : Hom(^) i—> Harm(£) is 
effectively a projection operator (which is not inversible). 
An equivalence relation 

We may define an equivalence relation in Hom(i'): , 

P^P'<^P-P' = r^ Q; QGHom(£- 2). (11) 

Two £ - homogeneous polynomials are equivalent if their difference is a 
multiple of the monomial = X ■ X. Then it is easy to check that 
two polynomials are equivalent iff they have the same harmonic projec- 
tion. Thus the vector space Harm(^) appears as the quotient of Hom(^) 
through the equivalence relation: each ^-harmonic polynomial represents 
an equivalence class in Hom(^). 

In tensorial notations, the harmonic projector 11 takes the traceless 
part of the tensor (see, e.g., 0). This allows to infer 

H = Il{vecH). (12) 
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All these relations are similarly verified by polynomials with real co- 
efficients. 

The correspondence being one to one (Sylvester's theorem), we may 
invert [the reduction to Vec(^, IR) of] the 11 operator. This provides the 
interpretation of CHS's result as the action of n~^: 

Harm(£, IR) ^ Vec(£, IR) 
n : H = UvecH ^ vecH, 
n"^ : H ^ vecH = n'^H. (13) 

This allows to see also Vec(^, IR) as the set of equivalence classes in 
Hom(^, IR). Note however that Vec(^, IR) is not a subvector space of 
Hom(^, IR). But we may provide a vector space structure to Vec(^, IR) 
by defining the " sum " : 

P + Q = n [U-^P + U-^Q] , P,Q(£ Vec(^,IR). (14) 

For instance, we have + y'^ = ~x^. 



2.4 Link with the multipole coefficients 

The CHS's formula allows to characterize any multipole /(f) by the con- 
stant A and the coordinates of the £ unit vectors. It is interesting to make 
the link with the other description given by the 21 + 1 coefficients ftm of 
the usual harmonic decomposition Q. 
We start from the identity 

ai a2 ...ae = ^ ^ {ei ai +€2 a2 + ... + ei a^Y ; B ^def ^! ^ 1, 

ei ei ...Ef £1 El ...Ef 

(15) 

where each takes the two values -1 and 1. It implies 



{ui- X) {u2- X) ...{ue ■ X) = ^ ^ [{ei ui +62 U2 + ... + ei ue) ■ XY . 

ei ei ...El 

(16) 

Assuming unit vectors, this may be rewritten 

(ei COS^l 4-62 COS02 + ... + Q cosOiY, 

ei ei ...eg 

with cos6'i — Ui ■ X. From the calculation in the appendix, it results that 
n[{uj ■ X) {U2 ■ X) ...{ue ■X)]=J2 "^"^ Yinri^), 

m 

with 

aim^^ ^ yftn(ei Ml + €2 U2 + ... + et Ui), (17) 
ei ei ...€f 

with Ki given by l|82|l . 

This formula gives the coefficients of the harmonic decomposition as a 
function of the multipole vectors. The reciprocal calculation can be made 
in the tensorial formalism as indicated by CHS. 
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3 Extension to homogeneous Polynomi- 
als 

By an arbitrary ^— homogeneous polynomial P G Hom(£, IR) is uniquely 
projected to HP e Harm(^, IR), with HP = P+r^ Q', Q' G Hom(^-2,]R). 
On the other hand, CHS result implies 

nP = vec{UP)+r^ Q", vec{nP) € Vec(£,IR), Q" £ Hom(£-2,IR). (18) 

These two relation imply 

P = n^^nP + r^ Q; Q e Hom(£- 2,IR). (19) 

This extension of CHS's result to Hom(^, IR) had been found by KW. 
The present derivation provides a shorter demonstration: any real £— 
homogeneous polynomial P G Hom(^, IR) has the unique decomposition 

P = vecP + Q, (20) 

with vecP G Vec(£, IR) and Q G Hom(^- 2, IR). The CHS's result appears 
as a special case. 

3.1 Stability 

For any function f on S^, like for instance the temperature of the CMB, 
let us call 

L 

f<L = Yl fw (21) 

1=0 

its approximation by the L first multipoles. For each value of L, f<L may 
be developed as a sum of vectorial polynomials, according to the Corrolary 
2 of KW. Do the vectorial polynomials in the development change when 
L increases ? This is the stability problem asked by CHS and KW (their 
section VI: " an open question "), to which we answer here. 

To answer the question, let us consider two successive approximations 
f<L and /<L+i of the same function /. According to the KW's Corol- 
lary 2, they can be developed uniquely as 

L L+1 

f<L =Y Ve and f<L+i = ^ We, 
e=o £=0 

where Ve, Wt G Vec(£). As quoted by KW, Wl+i + Vl+i = 0, Wl = Vl 
and Wl-1 Vl-i in general. The question of stability concerns the 
possible equality between the other Ve and We, for i < L — 1. 
Let us assume stability, i.e., that 

Ve = We,W<L~ 1. (22) 

This would imply 

/(i+i) = f<L+i - f<L = Wl+1 + Wl-1 - Vl-1. (23) 
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First, let us note that, for a given f<L fixed, f(L+i) can be [the reduc- 
tion to S'^ of] any (L + l)-harmonic polynomial. On the other hand, the 
relation 12311 would hold for any /<l, and in particular when Vl-i = 0. 
This would imply that = Wl+i + Wl-i would hold for any (L + 1)- 

harmonic polynomial. In other words, 

Vi/ e Harm(L + l), H = V1 + V2, Vi G Vec(L + l), V2 £ Vec(L-l) : (24) 

the decomposition, which is unique, would involve two terms only. This is 
clearly wrong and, thus, the stability hypothesis 12211 is not true: no term 
in the decomposition is stable (excepted in the very special case, where all 
the monomial on the decomposition of P are already in vectorial form). 
Expansion of the exponential 

To judge the severity of this instability, it may be convenient to exam- 
ine a simple example, namely the case of the exponential f{x) = exp{k.x), 
where k and x are vectors of IR"^. The usual development of the expo- 
nential gives in fact its exact multipole vector expansion (in this case, as 
expected, the unit vector k = k/\k\— k/K is the unique one appearing): 

/w = E ^ (^-^r- (25) 

^ — ' n! 

n 

On the other hand, the well known multipole decomposition of the 
exponential implies (on the sphere) 

/(^)(x) = {2e + l) / MK) Pti-ik-x). (26) 

Here, je is the spherical Bessel function, and the Legendre Polynomials 
Pe are extended to complex arguments, insuring that /(^) {x) takes real 
values. It results 

L 

f<L{x)= E (2^ + 1) / j,(Jf) Pf(-jfc-a;). (27) 

The usual expansion of the Legendre Polynomial allows to estimate the 
higher (L) order term in this development, as 

(2L + l)i^(7f) ^i^(fe.x)^ (28) 

(the imaginers disappear). 

This is the leading term in the decomposition of /<l according to the 
corollary of KW. Thus, in this simple case, we may compare the term of 
order L in the expansion of the " complete " function /, given by l|25|l . 
with the corresponding term 12811 in the order L approximation f<L- Their 
ratio is given by 

jlCFO (2L + 1)! 
~ A'i 2^ (L!) ■ 
It differs from 1, but it tends towards unity when £ goes to infinity, as 
shown in Figure Q (in the case K — 1). 
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3.2 A stable decomposition 



This instability motivates the search for an other decomposition, which 
is in fact provided by the MaxweU representation . Starting from the 
multipole expansion 12111 of a function /, we simply replace each element 
by its Maxwell representation. This gives the decomposition 

f<L{x) = I^Ao r + Ai r'^ Vui,i + ... + Af r'^^'^^ .i-'-V^^, + ... 

+ r'^+' ,,...V„^ (i), (29) 

which is stable by construction. (I thank Jeff Weeks for the idea of this 
concise demonstration) . 

3.3 Anisotropy 

The multipole vectors offer an immediate advantage: they are rotated by 
S0(3) as vectors, i.e., according to the usual vector representation. Thus, 
they seem very appropriate to study the possible anisotropy of the CMB 
data. 

For instance, when T is a Gaussian isotropic random field, all the 
Tf^) are independent functions. Thus no correlation should exist between 
vectors corresponding to different multipole values I. The situation is 
not so simple for the different vectors in the decomposition of a given 
multipole. In the absence of any preferred direction, no peculiar vector 
could appear in the decomposition. Thus, for any realization, the vectors 
obtained should be isotropically distributed. For the quadrupole, H| have 
studied the algebraic and statistical independences of the coefficients in 
this decomposition, for a Gaussian random process. 

This gives a high significance to the results that CHS and KW report, 
from their analyses, of an " astonishing " quadrupole and octopole align- 
ment. This clear breaking of isotropy sets the question to interpret it 
as a chance effect, a contamination of the data, or a cosmological effect. 
In this latter regard, it is tempting to invoke an universe with multi- 
connected spatial topology, where the large scale isotropy disappears, and 
is replaced by partial isotropy, i.e., symmetry under a group H related to 
the holonomy group F. It is however not so clear that this may explain 
these results. 

It is well known that a multi-connected space imposes a quantifica- 
tion of the wave vectors appearing in the mode decomposition of spatial 
fluctuations, which is reflected in the distribution of angular temperature 
fluctuations of the CMB. The statistical distribution of the latter must 
therefore be invariant under H. This has for consequence that, for a given 
multipole, the distribution of the vector Ui must be iJ-invariant, with the 
immediate consequence that the number of such vectors must be an in- 
teger divisor of the order uh of H . This implies that, for a perfectly 
representative distribution, the vectors Ui can be present only if I (or a 
divisor of t) divides uh- For instance, if the quadrupole and octopole 
alignment is significative in this regard, this would imply that uh is a 
multiple of 6, which puts strong constraints on H , and thus on F. 
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The question of anisotropy has also been treated by 0|. They pro- 
pose to extract from the set of vectors Ui (and the constant ^4) related 
to a given value of I, a set of 2^ — 2 rotation invariant quantities plus 3 
rotational degrees of freedom. This interesting suggestion, however, only 
provides a partial answer to the question of anisotropy. Firstly, the rota- 
tion invariant quantities provide information on anisotropy. For example, 
a multi-connected space will impose some fixed angles between the vec- 
tors, in order that the set of vectors be H invariant (for instance, right 
angles in a toroidal universe). Thus, even if they are rotation invariant, 
the scalar products between the vectors bring valuable information about 
anisotropy. On the other hand, the examination of (for instance) the ro- 
tation properties of " anchor vectors " when I increases, has no simple 
interpretation. For instance, in the simple case of a toroidal universe, the 
different multipoles (different ^'s) may favor vectors which are the sides, or 
the diagonals, or different preferred vectors. In this case, despite a strong 
and well defined anisotropy, the anchor vectors would show no alignment 
(although they would exhibit some specific and well defined correlations 
which, however, cannot only be predicted a priori with a specified partic- 
ular model). Thus we conclude that the interpretation of the multipole 
vectors in terms of isotropy/anisotropy of the data remains an open ques- 
tion. 

In fact (excepted for I = 1,2), there is no standard way to associate 
to the I multipole a definite spatial orientation (which belongs to the 3 
dimensional representation of SO (3)). The only exceptions are the dipole 
{(. = 1), which selects one preferred direction; and the quadrupole, which 
selects a preferred orthonormal frame (see also |2] for an interpretation 
of the dipole and quadrupole vectors). In fact, an I multipole may be 
seen as a function on the fuzzy (non commutative) sphere Sfuzzy,2e+i, 
which is itself an approximation of the ordinary sphere by a set of 21+ 1 
cells. 

4 Extension to complex polynomials 

The results of CHS and KW were obtained for polynomials, and vectors 
in the decomposition, with real coefficients. Since, for instance, the usual 
spherical harmonics correspond to polynomials with complex (not real) 
coefficients, it seems interesting to try to generalize the previous results. 
(Note that a basis of real spherical harmonics ofi'ers a limited interest, 
because they are not eigenvectors of the projected angular momentum). 

The KW s demonstration is based on the existence of 21 common roots 
(in CP^), to the two equations P{X) = and — 0. This remains true 
when the coefficients of P{X) are complex. However, in this case, the 
roots are no more two by two complex conjugates. Although a similar 
demonstration can be performed, there are two main difi'erences : 

- it remains true that the 21 roots can be grouped by pairs, to gen- 
erate i lines. But there is no canonical way to do it, as in the real case 
(associating the complex conjugates). The consequence is that the result- 
ing decomposition is not unique: there are A*" different decompositions, N 
being the number of distinct ways to group the 21 roots (which may be 
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distinct or not) in pairs. 

- In general, the roots are not complex conjugates. Thus the constructed 
lines are not real. This means that the vectors obtained in the decompo- 
sition have complex coefficients. 

We obtain thus the following result : 

An i!-homogcncous polynomial with complex cocfficionts, P G Hom(f) 
can be decomposed in a finite rmmbor of different ways, under the form 

P{X) {X ■vi)...{X ■ve)+X'' Q, (30) 

where the Vi have now complex coefficients and Q € Hom(^ — 2). Given 
the normalization constant a £ C, it is not a restriction to assume these 
vectors unitary = I) or null {v'^ — 0). 

When P has real coefficients, one decomposition (among the others) is 
canonical, and involves only real (unit) vectors : this is the KW's result. 
There are however, in general, other decomposition involving complex 
vectors, two by two conjugates. 

4.1 Some examples 

• P = has the real decomposition P = —z^ + r^, but also the 
complex one P = [x + iy) {x — iy) . 

• P = x^ — y^ has the real decomposition P = {x + y)(x — y) , but also 
the complex one P = —{2y'^+z'^)+r'^ = —{^/2 y—iz){\/2 y+iz)+r^ . 

• For the toy quadrupole of KW, their equation (19), there are, as 
expected, the 4 roots given by their equation (22). By associat- 
ing the complex conjugates, they obtained their decomposition (26) 
with real vectors. But different associations provide the two other 

decompositions: 

l/6[x{l-iV2)+Sy+z {l-iV2)] [x{l+iV2)+3 y+z {l+iV2)]+3r'^ /2 
and 

-~l/2[-x{l+iV2)+y+z {-l+iV2)] [x{-l+iV2)+y-z {l+iV2)]+r^ /2. 

• The toy octopole of KW (27) gives the 6 roots (29,30). The associ- 
ation of the complex conjugates gives the decomposition (32) with 
real vectors. The different groupings give 

{i x + y - z) {-i x + y - z) /2 + {y - z/2), 

and 

-{ix + y + z) {-ix + y + z) /2 + {y + z/2). 

4.2 The spherical hcirmonics 

Note that the operators Ji = —ietjk dk, and thus J^, commute with A 
and with 11. From this, it results that the spherical harmonics are given 

by 

Ye,^=Uz'-"' {x + iyr. (31) 



10 



In particular, 

Yu =n {x + iyf = {x + iy)\ 

Yu-x = n z (a; + iyf-^ ^ z (x + iyf''', 

Yu-2 = n «2 (x + iyf-^ = z' (x + iyy-^ - rVs {x + iyf 

y,o = n/ = 4^p,(.), 



5 Conclusion 

As proven by previous results, the multipole vector decomposition offers 
a new and promising technique. However, the instability of the decom- 
position, and the difficulty to interpret anisotropy, clearly demand deeper 
analysis. In this regard, the result presented here, namely the link with the 
Maxwell representation, the interpretation in terms of harmonic projec- 
tion, the established correspondence with the usual harmonic development 
give some new perspectives. 

It remains to analyze the practical relevance of the new stable decom- 
position proposed here; also, the link between anisotropy and multipole 
vectors deserves further exploration, in particular in the frame of cosmo- 
logical models with non trivial topology. 

Moreover, the scalar product X • m, as a function of X, shows some 
similarities with a wavelet on the sphere. Thus, the multipole vector 
decomposition shows analogies with a product of wavelets. Since there 
is presently a lot of interest toward wavelets on the sphere, a further 
exploration of these analogies may appear fruitful. The construction of a 
new formalism using the multipole vector decomposition appears therefor 
as a very promising task. 



6 Appendix 

By definition, n(« ■ Xf , is ^—harmonic and can thus be expanded in 
spherical harmonics as X)^a^m(^) yim{X), for u and X on the sphere. 
Using spherical coordinates, we have 

X = {cosO, sinO cosip, sin^ sinc^) and u = {cosO', sinO' cos (^', sin 0' sin<^'). 

Thus, n(X • uY = Tl[cos9 cos 9' + sinO sin 5' cos{(p — <p')] depends on <p' 
only through if — if' . Using Yim{X) = K(m. e""''' Pem{cosd), where the 
Kirn's are the usual normalization constants and the Pcm's the Legendre 
functions, this implies 

n{X.uY = ^ Ferr^iO') e-*""^' Y,m{X) = ^ e'"'*'^-'^'' Ke^ P(m{cose). 

m m 

By symmetry, wc have also 

n(X ■uf = J2 Fem{e) e^-(-^+'^') K,m Pem{cose'). 
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This implies, using symmetry arguments, 



n(X • uf = KeJ2 ■^«*'"(") YtrniX). 

m 

The normalization constant 



- {21+1) (2^!) (^2) 



can be calculated from Q, using the polynomial development of Pi 
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Figure 1: The ratio i? as a function of the muhipole index t, ior K = 1 
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